Single Photon Ignition of Two-photon Super-fluorescence 
through the Vacuum of Electromagnetic Field 

Nicolae A. Enaki 

q". Academiei str.5, Institute of Applied Physics of Academy of sciences of Moldova, 

O ' Chisinau, MD 2028, Republic of Moldova 

(N ■ 

O : December 8, 2010 

<D ■ 

Q: 

C*"** \ Abstract 

The ignition of two-quantum collective emission of inverted sub-ensemble of radiators due to 
mutual interaction of this sub-ensemble with other two dipole active atomic subsystems in process 
of two-photon exchanges between the atoms through the vacuum field is proposed. The three 

-f— > ■ 

£^ ■ particle resonances between two-photon and single quantum transitions of inverted radiators from 

a ■ 

r-{ • the ensemble are proposed for acceleration of collective decay rate of bi-photons, obtained relatively 

dipole- forbidden transitions of excited atomic sub-ensemble. This mutual interaction between three 
super-fluorescent processes in subatomic ensembles take place relatively dipole-forbidden transitions 
in one of radiator subsystem. The collective resonance emission and absorption of two-quanta 

IT) ' 

' have nontraditional behavior, accompanied with acceleration and inhibition of collective emission 

processes of photons. 

m ; 
(N ; 

1 Introduction 

T— I . 

A great deal of attention is currently devoted to the problem of coherence which appears not only between 
the quanta but between groups of quanta too. The generation of non-classical coherent electromagnetic 

b : 

field in multi-photon emission and the interaction of coherent radiation with matter (nuclei, atoms and 
solids) have been subjects of a number of theoretical and experimental studies in recent years PQ-[3]. 
Examples include the higher-order coherence in multi-photon generation of light the two-photon micro- 
maser emission [3] , two-photon lasers the parametric down conversion, four- wave mixing and other effects 
in optical diapason [2 , and the possibility of coherent generation of photons in x-ray and gamma- ray 
spectral regions. 

In this article it is proposed to investigate the cooperative two-photon emission from inverted system 
of radiators stimulated by single photon super-fluorescent pulses in two-quantum resonance with dipole 
forbidden atomic transition. Since the two-photon cooperative phenomenon has the small two-photon 
cooperative emission time [5] , we propose to extend our attention to the new tape of cooperative res- 
onance interaction between three radiators in which single photon transitions of two radiators which 
enter in two-photon resonance with dipole forbidden transition of third atom. This cooperative three 
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particle interaction take place through the vacuum fluctuations of electromagnetic field and can amplify 
or diminish the spontaneous emission rates of the atoms. In order to obtain more powerful pulses of 
entangled photons it is proposed the cooperative interaction between three atomic subsystems in which 
one of them are inverted relatively dipole forbidden transition |2S" > — > of Hydrogen like or He- 
lium Like atoms [l]-[8]. Taking in to account the elementary acts of two photon interaction between 
radiators we archived the improvement of two-photon emission rate of the system of radiators in com- 
parison with two-photon super- fluorescence [5] . In this article it is examined the mutual influence of two 
single-photon super-fluorescence processes and two-photon cooperative emission of the atomic system 
relatively dipole forbidden transition. The phenomenon of new cooperative emission takes in to account 
the three particle mutual interaction with vacuum of electromagnetic field in which the product of vac- 
uum polarization of two atoms enter in to resonance with two-photon polarization of dipole forbidden 
transition of Hydrogen- like or Helium-like radiator. It has been shown that in the process of spontaneous 
radiation, the radiators (nuclei, atoms) enter a regime of single and two-photon super-radiance and the 
rate of photon pair (bi-photon) emission increases (or decreases) due to new three particle cooperative 
phenomenon, which appear between single and two-photon spontaneous emission subgroups of radia- 
tors. It has been demonstrated, that for hydrogen-like and helium-like atoms [10] the dipole-forbidden 
transitions can generate more powerful pulse of entangled photon pairs (bi-photons) under the influence 
of single photon super-radiance. 

It is important to note that for coherent radiation of such system, was studied for the dimension 
of a radiating system smaller than the radiation wavelength. It is, however, interesting to study this 
type of cooperative emission between three radiator subsystems in extended system of radiators. The 
possibilities of two-photon cooperative resonances between three radiators replaced at distance larger 
than emission wavelength are studied too. I emphasize here, that the problem of cooperation between 
two single photon cooperative emission subsystems and one two-photon cooperative emission subsystem 
is more complicated than the similar problem of single 11 lj or two-photon if)] super-radiances in extended 
system. In Dicke's super-radiance, the exchange integral between j-th and l-th atoms is described by 
more simple exchange integral proportional to sin[kQrji]/(korji) while in two-photon super-radiance by 
more complicated function sin[(2/co — k)rji]/[(2ko — k)rji] x aisx[krji]/(krji), where Tji and ki are the 
distance between the radiators and wave vector of emitted photons respectively. The three particle 
exchange integral between j-th, m-th and Z-th atoms was obtained in this paper taking in to account 
the two-quantum exchanges between two radiators proposed in papers [9], |12j . The more complicated 
exchange integrals between two radiators with dipole active transition and one radiator with dipole 
forbidden transition is given in Appendix of this paper. 

2 Interaction Hamiltonian and Master Equation 

Let us consider the interaction of three subsystems of radiators R, S, and D thorough vacuum of 
electromagnetic field. The first two groups, R and S, are prepared in excited state \e r ) ® \e s ) and can 
pass in to Decke super-radiance regime relatively the dipole active transitions e r — > g r and e s — > g s 
at frequencies u> r and uj s (see figure [I}. The D atomic subsystem is prepared in excited state |e<j) and 



2 



R 




2<>>o =co, +w s 



Figure 1: The resonance between two-photon transitions of D atomic subsystem and two dipole active 
atomic subsystemsi? and S . As an example is represented three atoms D, R and S situated at relatively 
distances r<j s , r ( i r and r rs . One of condition of exchange energies between the subsystems is the e 
resonance between two-photon and single photon transitions 2loq = uj r + lu s . 

relatively dipole forbidden transition a — > go, and can pass in the ground state \ga) simultaneously 
generation two quanta [5]. Let us consider the simple cooperative stimulation of two-photon emission 
of D system stimulated by R and S radiator subsystems. 

In this case it is established the resonance between the dipole active atomic subgroups R, S and 
dipole forbidden radiators of D ensemble and the cooperative stimulation of two-quantum collective 
transition is possible. Indeed, considering that the conservation energy law is established between these 
groups, h(u) r +u s ) = 2hwQ, one can proposed the following Hamiltonian of interaction of radiators with 
electromagnetic field 

H = Hq + XHj; 

N r N s N 

H = hhj k a\a k + ftw r Rzj + ^ t^sSzi + 2 hu> D zm ; 

k 3 = 1 1=1 rn=l 

N a 

A-fff = -y^y^(d r ,g fc ){i?+a fe exp[i(k,r j )] +i?^4exp[-i(k,r j )]} 
k 3 = 1 
N b 

- y^y^(d s ,g fe ){S' ; + a fc exp[i(k,r ; )] + Sfa\ exp[-i(k, rj)]} 

k 1=1 
N b 

- ^2 y^(n e g,e Al )(n ej ;,e A2 )(7(a;i,W2) 

k 1: k 2 m=l 

x {D+ n a k2 a kl exp[i(ki +k 2 ,r m )] 

,t „t 



An4i4 2 ex PH( k i+ k 2> r 0]}- (1) 



Here 



, \ ^23^315/^2^/0! f 1 1 | \2ltfuX>k 

2h y uj 32 + Wfej W31 - ujk 2 J v v 

afc and at are annihilation and creation operators of EMF photons with wave vector k, polarization e x 
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and frequency uik', d r and d s are dipole momentum transition between the ground and excited states 
for R and S atomic subsystems; d e i and d eg are dipole momentum transitions in the three level system 
of atomic group D. The operators of i?, S, and D atomic subsystems satisfy the commutation relations 
for SU(2) algebra [J+, J~] = 2J Z ; [J, , J±] = ±J ± , where J± is equivalent with R ± , S± and D ± . 
Invertin operator J z is consider similar to R z , S z and D z respectively. The operators of electromagnetic 
field satisfy the commutation relation [a^a^,,] = 5 k k i ; [aj.,a^,,] = 0, where fc = (k,A) is the wave vector 
and polarization of the photon. Taking in to account the Hamiltonian [TJ let us represent the solution 
of Haisenberg equation through the sources and free part operators 

dk(i) = a k {0) exp[-iw fc t] + a ks (t), (2) 

where the source part is 

a ks (t) = y-^-^exp[-i(k, r ; ) / dr exp[-iaj k T]R~ (t - r) 

+i exp[-i(k, Tj) / dr exp[-iu> k T]S~ (t - r) 



3=1 



+2 . ^ ^ k n efl ,e Al )(n ei ,e A ) g (a, 1 ,a,) exp[ _ j(ki + 

n—l 



I 



dTexp[-iu k T]D n (t - T)a ] ki {t - r); [a\ s {t) = [a s (t)X 
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Taking in to account that Ofc(0)|0) p / l = (0| p /ial(0) = we can partially eliminate the EMF field operators 
from the mean value of Hesenberg equation for arbitrary atomic operator 0(t) 

|(0(i)} = -i^^^pl([ J R;(i),0(«)]a fcs (t))exp[i(k,r J -)] 

k j=l 1 



E E ^T^l^W. 0(*)]a*. (*)) exp[i(k, r,)] 



h 

k 1=1 

Eyv (n eg , e Al (fci))(n ei , e A (fc))g(cj fcl , ajfc) 
^ h 

ki,k m—1 

x ([£>+(«), 0(t)}a kl (t)a ks (t)) exp^ + k, r m )] + H.C.(0+ -► O). (3) 

Here the mean values of Hesenberg operators are considered taking into account the initial state of 
the system |5' r (0)) <E) \0}phr where |\I/ r (0)) is the state of radiator subsystem, and \0} p h is the vacuum 
state of EMF. We are interested in the total elimination of operators of electromagnetic field from the 
expression ([3]) For elimination of operators of electromagnetic field we formulate the lemma 

Lemma 1 If Bose a k (t) and a k (t) operators lie between the two operators of the atomic subsystem 
A(ti) and Bfo) (A(t{) , -6(^2) don't contain the operators a k and a k ) belonging to other times, the 
elimination of the free part of these operators yields the following expression for the correlation: 

{A{tx)a k {t)B{t 2 )) = (A(ti)o fc ,(t)B(t 2 )> 
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e- to *(*-fc>(A(i 1 )[a fca (t 2 ),B(t 2 )]), 

(A( tl )a+(t)B(t 2 )) = (A{ tl )a + ks {t)(t)B{t2)) 

- e^~^ {[A^atMmt^)) ■ (4) 

Proof. The commutations in Q play the highest role in the two-photon spontaneous emission and only 
such commutations bring the main contribution to the two-photon process. The problem is reduced to 
the elimination of vacuum part lies between the operators A(ti) and B(i,2) 

< A{tx)a k {t)B{t2) >=< A{t x ){al{t) + a ks {t))B{t 2 ) > (5) 

Since a k (t) = a k (t) + a ks (t), we will represent the vacuum part a k (t) = a k (0) exp[— iui k t] through the 
vacuum-operator at time t\ and taking in to account the identity @ we can represent the vacuum part 
in the following form a%(t) = a£(t 2 )e _;iWfc(t ~* 2) = {a k (t 2 ) - a ks {t 2 )}e,- luJk(t - t2 \ After substitution of 
a k (t) into the correlation it is obtain 

{A{tx)a k {t)B{t 2 )) = {A{tx)a ks {t)B{t 2 )) 

+ e-^-^iAih^akih) - a ks (t 2 )}B(t 2 )), 

We observe that <Xfc(<2) commutes with the operator B(t 2 ). Consequently taking into account that 
cik(t 2 ) B(t 2 )\0 >= B(t 2 )a ks (t 2 )\0 >, it is easily obtain that 
(A(ii){a k (t2) - a ks (t 2 )} B{t 2 )) = - (B(t 2 )[a ks (t 2 ), B(t 2 )}). 
This relation proofs the Lemma. ■ 

This lemma Q can be used in the last term of generalized equation ([3]) for correlation functions 
([Df(t),0(t)}a k (t)S-(t-T)) , ([D+(t),0(t)]a k (t)Ri(t-r)) and ([Rf(t),0(t)}D-(t-r)al(t-r)}, 
([Sf(t),0(t)]D~(t — T)a ki (t — r)). Indeed taking in to account the lemma ((4]) the above correlation 
functions can be represented through atomic operators 

([D+(t),0(t)]a k (t)S-(t-r)) = ([Df(t),0(t)]a ks (t)S;(t - r)> 

e-^([D+(t),0(t)][a ks (t-r),S-(t-r)]), (6) 

([R+(t),0(t)}al(t -T)D-(t-r)) = ([R+ (t),0{t)]a\ s {t - r)D~{t - r)) 

e-^({{R+(t),0(t)},ai(t)}Dn(t-T)). (7) 

The interaction between the atomic subsystems can be found in the third order of interaction constants 
with the subsystems S, R and D respectively (d r , g k )(d s , g k )q(uji, uj 2 ) According with this condition 
the smooth correlation functions is obtained only for the following terms of expressions and ([7]) : 
(R[{t - T')[D+(t),0(t)]S- (t - t)) and (R+(t - T')[R+(t),0(t)]D-(t - t)). The contribution of other 
terms of the expressions ([B]) and ([7]) give the contribution more hair order on the decomposition on the 
small parameter A of the interaction Hamiltonian ([TJ . 

The lemma Q is non-applicable for correlation functions in which it is meet simultaneously the cre- 
ation and annihilations Boson operators belonging to different time intervals: (\Dj~(t), 0(i)]a k (i)D~ (t — 
T ) a k 1 T )) an d its hermit conjugate part (a kl (t — r) D+ (t — T)a\(t)[0(t) , DJ (t)]). In order to eliminate 
the vacuum part of operators a k (t) and a k (t') let us formulate the following rule. 
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Lemma 2 If the operators A(t±) contains the creation operators of EMF and B{t 2 ) contains the annihi- 
lation operators of EMF the elimination of vacuum part of annihilation ak (t) or creation at (t) operators 
situated between these operators A{t\) and B(t 2 ) takes place according with Lemma 1. In opposite case, 
when operator A(t\) can be represented through the product of atomic operator A{t\) and annihilation 
field operators A(ti)=A(ti)a kl (ti)a k2 (ti)...a kn (ti) . the operator B{t 2 ) is represented through the prod- 
uct of creation field operators and atomic operator B(t 2 ) so that B(t 2 )=B(t 2 )a ki (t 2 )a k2 (t 2 )...a k (£2) the 
elimination of vacuum part of the operators a k (t) and at (t) can be represented in the following form 

(A(jti)a k (jt)B{t2)) = (A(tx)a ka (t)B(t 2 )) 
- exp[-iu k (t - t 2 )]{<A(ti)[o fca (t 2 ), B{t 2 )}) 
-5 kM (A(t 1 )B(t 2 )al(t 2 )...aijt 2 )) 

- h,k m (A(ti)B(t 2 )a+ 1 (i a )ot i (t a )... a t m _ 1 (i a ))} J (8) 

(A( tl )al(t)B(t 2 )) = (A( tl )al(t)B(t 2 )) 
-exp[ioj k (t-t 2 )]{([A(ti)al s (ti)}B(t 2 )} 
S k , kl (A(t 1 )a k2 (t 1 )...a kn {t 1 )B(t 2 )) 

- 5 k , kn (A(h)a kl (ii)a fe2 (ti)...^., (ti)B(i 2 )». (9) 

Proof. Taking in to account the lemma (@|, we can represent the third correlation (A(ti)ak(t)B(t 2 )) 
of expression ((5J in the following form 

{A{tx)a k {t)B{t 2 )) = {A{t x )a ks {t)B{t 2 )) 

+ exp[-iu) k (t - t 2 )] (A(ti)(a k (t 2 ) - a ks (t 2 ))B(t 2 )}) 

- eBtp[-MJ fc (t-t2)]{<A(t 1 )[a k ,(t 2 ),B(*2)]). (10) 

According with explicit expression of operator, B(t 2 ) — B{t 2 )a\. {t 2 )a' k (t 2 )...a\, (t 2 ), let us introduced 
it in the third term of right hand site of expression (|TU)) . Following the commutation roles of boson 
operators of electromagnetic field, the operator a k (t 2 ) can be permuted in the right hand site of the 
correlation. Taking in to consideration that (a ks (t 2 ) + a kv (t 2 )) |0) = ak s (t 2 ) |0) , this term becomes 

(A(ti)o fc (ia)B(t 2 )oi i (t2)ai a (ia)...oi ra (ta)) = 

5 kM (A(t 1 )B(t 2 )a{ 2 (t 2 )...aijt 2 )) 

+■■■ + h,k m (A(i 1 )6(i 2 )4 i (t 2 )al (t 2 )...al m i (£ 2 )) 

+ (^(t 1 )B(i 2 )4 i (£2)4 2 (^)-4 m (t2)a fes (t2)). (11) 

Introducing this relation in (jlip it is not difficult to observe that the new expression for correlation 
(A(ti)ak(t)B(t 2 )) coincides with ([8]). The similar procedure of permutation of vacuum part of creation 
operator a\(t) demonstrates the identity (j9]) of Lemma 2. ■ 

According with lemma © we obtain the following expression for correlation function 
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([D+(t),0(t)}a k (t)D-(t - T )al(t - r)> = ([D+(t),0(t)]a ks (t)D-(t - r^t - r)) 
-exp[-iu; k r\{([D+(t),0(t)}[a ks (t-T),D-(t - rja^t-r)]) - 5 fclfcl ([Z?+(i),0(i)]Z?-(t - r))()L2) 

The next step of elimination of operator a\. (t — r) from this expression must be do taking in-to account 
the lemma (j4]). When the first and second order interaction constants have the same small magnitude 
A ~ (d s , gfc) ~ q{uJk t i Wfe), in Born approximation we take in to account only the last term of expression 
(|12p . the interference contribution of which is proportional to A 3 . As follows from the representation 
([2|) and ((3]) the procedure of elimination mast continue. Indeed introducing again this equation in the 
right hand cite of equation (|14[) we obtain the following master equation for arbitrary operator 0(t) in 
thread approximation on the interaction constant A 

d(0(t)) ^ ^ (d a , gk ) 2 



dt 



EE i -^W tL J dTcM-^ k T + l (k,r 3 -r l )]([R!(t),0(t)}Ri(t-T)) 

k l,j=l q 

E E / drexp[-iu k T + i(k, r, *)]([£$ (t), 0(t)]Sf(t r)> 

J„ 11 1 « 



k 1,1=1 

^ ^ (n e3 ,e A ^ (n ^e A jV(^ 2 ) j {t)]D T {t - r)) 

ki,k 2 1,3=1 







X 



x exp[-i(2w - wjfej - w fe2 )T] exp[i(ki + k 2 ,rj - r ; )] 

■ ( n e 9 ,e A (k))(n e4 ,e A (k'))(?(wfc,w fe O f /" 
+ 1 E E E E ^3 dr dr exp[i(k, r„ - rj) - zw fc r] 

fc,fc' n=l 3=1 Z=l g q 

exp[z(k', r„ - r,) - iw fc /T / ][(d„g fc ,)(d r ,g fc ) <[£>+(*), 0(t)]i?7( f - T )Sf(t - r')) 
+ (d r ,g k ,)(d s , Sk ){[D+(t),0(t)]Sr(t - r)R-{t - r')>] 

N t t 

■ ST Y^ Y^ Y^ (n eg ,eA(k))(n ei ,eA(k'))<z(wfc,WfcO /" , /" 

- ? EEEE is / dT dTexp[iu k T-i(k,r m -n)] 

k,k> m=l 3=1 (=1 o o 

x expHwfcr' - iu^r' + i(k', - r m )][(d fl , g*)(d r , g k >)(Sf{t - r)[Rf(t), 0(t)}D m (t - t')> 

+ (d s , Ek )(d r , Ek ,)(R+(t-r)[S+(t),0(t)]D-(t-r'))]+H.c.(0+ -+ O). (13) 

The traditional Born-Marcov approximation in the right hand site of equation ( 1131) give us the 
divergent functions. In order to understood this we approximate the right hand site of equation (|2|) with 
following expression 



(d g fc ) Na 
1=1 



=i 

(n es ,e Al )(n ei ,e A )q(w 1 ,cj) 



3=1 

Aft 
n=l ki 



xD-{t)ai i (t)cxp[-i(k 1 +k,r n )C{u k +uj kl -2ujo), (14) 
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in the Born-Marcovian approximation [15] , |13j . The small parameter in this approximation is the 
ratio of retardation time to cooperative spontaneous emission times of the subsystem, t/t% « 1. Here 
iQ(x) — iP/x + 7t8(x) is the Heitler function [13]. [14] . represents k-summation in analogy with Cauchy 
principal value |15) . Introducing the operators (I14[) in equation @ and eliminating the boson operators 
of EMF, it is obtain the following equation for operator 0(t) in Born-Marcov approximation 



k l,j=l 

x exp[z(k, r 7 - - ri)]i^*(cj r — <^k)+/l/ J 



(d a ,gfc)(d b ,g fc ) 



ft 2 

k lj=l 



x ([5t(t),0(t)]5 i -(t))exp[i(k,r J . -r,)]iC*(w.-w fc ) 
+ y^ y^ (n eg , e A ) 2 (n e » , ey ) 2 g 2 {u>k , Uk> ) 



h 2 

k,k' 1,3=1 



x ([D+{t),0(t)}D-(t)) exp[i(k- k', rj - - rj)]»C*(w« - - w*) 
+ ^ y-' y> y^ y> (d a , gfc' ) (d b , g fc ) (n eg , e A ) (n ei , ey )g(^fc , ^ ) 

fc,fc' n=l J=l i=l ^ 
x ([D+(t),0(t)}Rj(t)Sf(t)) 

x exp[i(k, r„ - r ; ) + i(k',r„ - r ; )]iC*(w r - w fc )<*(w s - 

_ ^y-^ ^ ^2^2 ^ da,Skl ^ db,Sk2 ^ ne9,eXl ^ nel,ex2 ^ Wl, ^ 2S) 

k,k' m—1 l—l j — 1 

x [(5+(t)[i?/(t),O(t)]£>-(t))iC*K-a; fe )iC(a; s -a; fc 
+ (Rt(t)[Sf(t), 0(t)}D- (t))iC (w. - w fc )iC(w r - Wfc/)] 

exp[i(k,r ) + i(k',r,) -i(k + k',r m )] + H.C.(0+ -> O). (15) 



In the right hand part of the equation (|15j) the third order terms contain the resonances between 
the single photon radiators A ,B and two-photon radiator D described by the correlation functions 
(S+(t)[R+(t),0(t)]D-(t)), ([D+(t),0(t)]Rj(t)Sf(t)) and (R+(t)[S+(t),0(t)]D m (t)). As it is observed 
from equation (|15[) . these terms contain the product of the functions [P/ (w b — Wfc)] [P/(uJb — u)y )] which 
describe the principal value in the integration procedure on the variables k and A''. It is not difficult 
to observe that these integrals become divergent expressions In order to avoid these divergence in 
Appendixl it is proposed the integration procedure which takes in to account the retardation between 
the radiators in the representation of right hand site of equation (j 13[) . According with Appendixl the 
right hand site of master equation (|13|) takes the following non-divergent form 



1,1=0 



!<°(*)> ^ 27 E xr(j,i)([Rt(t),o(t)]Rr(t)) + -L £ Xs (jj)([s+(t),o(t)}s^t)) 

r 1,3=1 

1 N 



1,3=1 

N N r N s 

+ t— E E E u & *. m)([D+(t),o(t)]Rj(t)sr(t)) 

Tbsr m=l {=1 1=0 
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N N r N. 

I 



E E E V & Lm)[{S+(t)[Rf(t), 0(t)\D-(t)) + (R+(t)[S+(t),0(t)}D m (t))} 



4l ~ sbr TO =l j=i 1=0 

+ H.C.{0 + ->• O). (16) 

Here the spontaneous emission r, and exchange integral between radiators j and I , Xi(i>0 are 
defined in expressions (|28|). ((32]), (|34| and (|37| ol Appendix. This equation can be used for description 
of interaction between the dipole forbidden and dipole active systems of radiators. 

3 Kinetic equations for correlation functions 

In order to found the correlation process between dipole forbidden transitions of D subsystem and dipole 
active transitions in S and R subsystems of radiators let us found the equations for arbitrary atomic 
correlation functions of subsystems of radiators. According with the generalized equation (II 6[) it is 
obtain the following chain of equation for atomic correlations 



L = l 

N N s 



^2[V(j,l.m)(S+(t)Rf(t)D-(t)) - V*(j,l.m)(D+(t)RT(t)Sf(t))},(17) 



4Tsbr m-i l-n 



rl 1 Ns 

Tf (s z i(t)) = -—Y,xs(i,p)[(s+(t)s;(t)) + x ;(i,p)(s;(t)s^(t))} 



dt x v " 2r s 



N N s 

+ T~ E Y^[V(j,Lm)(S+(t)Rj(t)D-(t)) V*(j,l.m)(D+(t)RJ(t)Sr(t))l (18) 



ATsbr m I / -u 



, 1 N 

-(D zn {t)) = -—J2lI>(3J)(Di(t)Di(t))+i;(j,l)(D+(t)D-(t)) 



dt x *" x " 2r b 

i=i 

N T N 

I 



Y,Yp{j,l,n){D+{t)R-j{t)S^(t)) - U*(j,l,n){S+(t)R+(t)D-(t)), (19) 



2Tsrb =1 1=0 



d 1 Nr 

-(R+(t)R-(t)) = -J2[Xr(l,n){Rf(t)R zl (t)R-(t))+ Xr (nJ){Rtm z Amr(t))] 



dt" 3 w 1 w/ r, 

n—1 



~ 5T E Y,[V(j,k,m){S+ (t)Rf(t)R zl (t)D-(t)) 

^'sbr , , „ 
m— 1 k— 

+ V*(j,k.m)(D+(t)R zl (t)RJ(t)S^(t))}; (20) 



, 1 N " 

j t (Sf(t)Sf(t)) = -J2lXr(l,n)(Sf(t)S zl (t)S-(t)) + Xr (nJ)(S+(t)S zj (t)Sf(t))] 



n=l 

N N r 



^^[y(fc,Lm)(i?+(t)5t(t)^(i)£>-(t)) 



2Tsbr m 1 /,-! 



V^i.mXC+W^^rWiirCt))]; (21) 
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1 N 

(D+(t)D-(t)) = - ^[/,a,m)< J D+(t) J D rf (t)£>-(t)}+7?(n ) m)< J D+(t)£> zn (t)£>r(t)>] 



dt n 

m—1 



N r N s 



EE^' k,m)(D+(t)D zl (t)RJ(t)Sk(t)) 

'bsr - i , n 

CT (i, fc, m) (St (t)Rf (t)D zn (t)D+ (<))]. (22) 



, 1 ^ 7Vs 

-i{D+(t)ST (t)Rk (*)> = — E E (■?'> m )< 5 « mt(t)D zm (t)Sr (W (*)) 



1 N 



n=l j=l 

V(kJ.n){D+(t)D-(t)R zk (t)Sf(t)Sr(t))\ 

+ * E [-*■ (i, (At (t)sr (*)) 

+ - Xr ( ]l k)(D+ l (t)R zk (t)Si{t)R-{t)) 

+ -It{m,l)(DUt)D zm (t)S-(t)R k -(t))] (23) 

Let us consider the interaction three different atoms in interaction through vacuum EMF. Introducing 
the exited numbers for the atomic subsystems (N a ) — (J Z a(t)) +0.5 (here J ■<-» S, R, D a = s, r, d ) 
and correlation function between the atoms (F) = i[(D + (t)S~ (t)R~ (*)) - (S+(t)R+(t)D- (*))] we can 
obtain the closed system of equations from the chain of equations (|17M23j) . Indeed considering that the 
distance between the radiators is smaller then radiation wavelength: 5ft {U(j, k,m)} — $t{V(j, k, m)} = 1 
, we obtain the following closed system of equation 



(XI Td T s b r 



— [6<iV s iV r iV d ) - 2(7V s 7V r ) - (7V s 7V d ) - (iV r iV d )], 



^-(N s N T N d ) = -[L + ± + L](N a N r N d ), 

dt T r T b T s 

^-(N s N r ) = -[— + -}{N s N r ), 

dt T r T s 

±(N s N d ) = -[± + !]( Na N b ), 

dt T r Tb 

i(N.N d ) = -[- + -}(N s N b ), 

dt T r Tf, 

in which the new correlation functions between the atomic excitation is introduced (N s N r N d ), (N s N r ), 
(N s N d ), and (N r N d ). This system of equation is exactly solvable. The solution is 
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1.5 



Figure 2: The decay law of inversion, D z (t), of dipole forbidden radiator stimulated by single photon 
decay processes of two radiators for following parameters of the system: relative decay rates of S- and 
R- atoms ti,/t s — Tb/T r = 5. 



(N s N r N d ) = exp(-At); (N s N r ) = exp(-Bt),(N s N d ) = exp(-Ct); (N d N r ) = exp(-Dt) 
(NiNj) = exp[-(- + — )t], i,j->a,r,b\ 

Ti Tj 

exp(-At/2) r 1 - exp[-At/2] 1- exp[-(B - T^)t/2] 

[ {)> ~ r bsr 1 A (B-t- 1 ) 

2(1 - cxp[-(C - r- 1 )t/2]) 2(1 - cxp[-(C - r^)t/2}) 

(C-Tr- 1 ) (C-r^ 1 ) J ' 
t 

(Ni) = exp[-t/Ti] ~ — ! exp[-(t - t')/ n ]{F{t')), i=s, r, d, 

7~bsr J 


where the collective rates are defined A = l/r r + 1/u + l/r s ; £? = l/r r + l/r s ; C = l/r& + l/r s ; 
D = 1/t;, + 1 /r s . The solution of this system of equation is plotted in figure [2] It is observed the 
influence of single photon transition on the two quanta transitions. This influence drastically depends 
on the cooperative rate 1/Tb sr - 

As follows from figure [3] the cooperative exchanges between the radiator accelerate the two-photon 
decay processes so that the the derivation —^f- achieved the maximal value decreasing after that till 
zero value. 

Neglecting the quantum fluctuations of inversion operators (R Zj ), (S z l)i (D zn ), dipole-dipole corre- 
lations between the same radiators (t)Rj (t)) , (Sf(t)S~(t)) and (D+(t)D~(t)) and between sub- 
systems (Sf (t)Rj (£)-D„ (t)), (D+(t).RJ {t)Sf (i)) we can de-correlated the chain of equations (|17H23j) in 
order to obtain the closed system of equations 

for J^S, R, D; (J+ (t)J zl (t)J m (t)) = (J zl (t))(J+ (t)J m (t)) j^l^m; 

(S+(t)R+(t)R zl (t)D-(t)) = (RamstWRfWD-Ww - S ld (S+(t)R+(t)D-(t))-, 
(D+(t)D zl (t)R-(t)S^(t)) = {D zl {t))(D+{t)RJ(t)S^{t))^ n - 5 Un {D+{t)RJ(t)S^{t))- 
(S+ (t)R+ (t)D zn (t)R~ (t)S^ (t)) = (D m (i)) (Rf (t)Rp (*)) (^ + (t)Sk (*)) ¥* W Vp 
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Figure 3: The two-photon decay rate dD z (t)/dt stimulated by single photon processes for same param- 
eter of the system as in figure O 

+5 j , p Si >k (D zn (t)(R zj (t) + 0.5)(S zl (t) + 0.5)); 
(D+ (t)D~ (t)S zl (t)R+ (t)R- (i)) = (S zl (t))(D+ (t)D~ (t)) (R+ (t)R- (t)) m ^ k 

5j, k {{D zn {t) + 0.5)S zl (t)(R ZJ (t) + 0.5)); 
(D+ (t)D~ (t)R ZJ (t)S+ (t)S^ (t)> = (R zj (*)) (D+ (t)D- (t)) (S+ (t)S* (t)) m ^ k 

5 m ,nSiA(D zn (t) + 0.5)R zj {t)(S zl (t) + 0.5)). (24) 

Taking in to account the de-correlation ([24]) , we obtain the following closed system of equations 
iMt) = --{Nr(N r + 2)/4-R 2 z +R z }--^—F, 
is,® = -^{N s (N s +2)/A-S 2 + S z }--^-F, 
±D z (t) = -L{N(N + 2)/A-D 2 z +D z }-—F; 

OX T r 7~sbr 

T+ F = [^[S z (t)-l] + ^[R z (t)-l] + ^[D z (t)-l]]F 
+ — [2D Z {N 2 /A - S 2 }{N 2 /A - R 2 } 

7~bsr 

+ R Z {N 2 /A - S 2 }{N 2 /A - D 2 Z } + S Z {N 2 /A - R 2 Z }{N 2 /A - D 2 } 
+ (AN s N r Nexp[-A * t] - N s * N r * exp[-Bt] - 0, 5N S N exp[-Ci] 
- 0.5N r Nexp[-Dt]]. 

From this system of equations follows the oscillatory behavior of the decay rate of the inversion 
D z . Taking in to account the following relative expressions of the decay rates we obtain the numerical 
simulation of the inversion D z and its derivative (see figure U and figure [5]). As follows from this 
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Figure 4: The time dependence of inversion D z (t) for dipole forbidden radiator subsystem for following 
values of the parameter of the systemmuber of atoms in the subsystems S, R and D are, N s — N r — 
N = 50 respectively; the relative decay times of the subsystems are 17, /t s = Tb/r r = 6 ; the coupling 
parameter of these three system is changed between and 1. The oscillatory behavior of decay rate is 
observed. 

system of equation the increasing of decay rate of two-photon spontaneous emission is possible under 
the influence of single photon cooperative emission of two atomic subsystems. In figures 2] is plotted 
the time dependence of the inversion < D z (t) > of dipole forbidden radiators as function of the relative 
coupled parameter between the radiators Tb/r sr b. 

The same dependence is represented in figure [5] for the intensity of two-photon emission proportional 
to —d < D z (t) > jdt. As follows from these plots it is observed the mutual influences between single and 
two-photon super-radiance processes of three particle interaction. This effect plays an important role in 
the collective decay process of the systems of radiators with the dimension smaller than wavelength. 

4 Conclusion 

In this paper the effective interaction between three radiator subsystems in two-photon resonance is 
found using the method of elimination of operators of vacuum field. The new cooperative interaction 
between dipole-forbidden atomic subsystem and two-dipole active subsystems of radiators was proposed. 
The master equation [TH1 which describes the energy dissipation from the system due to mutual interac- 
tion between the radiators through the vacuum of electromagnetic field, was obtained. Using the chain 
of equatior il7H231 which describes the cooperative interaction between three radiator subsystems, it is 
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Figure 5: The cooperative decay rate of inversion ^fi^ of dipole forbidden radiator subsystem for 
same values of the parameters of the system as in figure 21 The increasing of decay rate of bi-photons 
is observed. 



obtained the closed system of equations for three radiators. Neglecting the quantum fluctuation of the 
inversion, the de-correlation method of the this chain of equation is proposed [24] in order to describe 
numerically the behavior of mutual influences of single and two-photon super-radiance processes. As a 
consequence of effective interaction between three radiators through two-photon resonance processes of 
inverted systems increase substantially in process cooperative decay of the system The three particle 
exchange integral has been established and the influence of this effect on the behavior cooperative decay 
of the atomic subsystems was estimated (see figure 0] and figure [SJ . Similar experimental situation can 
bi realized in exited atomic (for example transitions in Cs atoms [17]) or nuclei (for example 193m jY, 
195m Pt and 103m Rh nuclei [16]) subsystems in resonance interaction through vacuum field. 



5 Appendix: Exchange integrals 

In order to estimate all exchange integrals in equation (|13p let us firstly found the well known exchange 
integral between two radiators in single photon interaction with vacuum of electromagnetic field. In the 
first terms of equation (|13|) the retardation can be found integration firstly on the k vector respectively 

ir OO t 

V-i = ( 27r )2 fe 3 J ^kduk J dQ k J dr exp[i(cJi - w*)t] 

1,3=1 o ' 

x(l- (ek,nd)){[Jj~(t),0(t)] Jf (t-r)) exp[iwfe rji cos 9], 
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"5 ? 



b. (25) 



Here the frequency ui^ corresponds to A and B atomic systems; i — r,s; for i — r operators Jj, J[~ 
corresponds to Rj , R~[ and for i = s these operators corresponds to Sj , Sf Passing to new variable 
v = LJk — <jJi and considering that the smooth function cok under integral can be approximation with 
uf, we obtain the following approximate expression of thirst order exchange integrals 

oo oo 

ujf J dv exp{iv (t — rji cos 9/ c)} ~ uif J dv exp[iv(r — rji cos 6/c)] 

—u>k — oo 

= 2nujfS(T - rji cos 9/ c). (26) 

uii is the emission frequency relatively the dipole active transitions of the R and S atomic subsystems. 
In this approximation I obtain the following integral on angle 9 and retardation r 

TT t 

V 3 \ = *M [d9 [ dTsm95(T - r ]l cos9/c)D l (9)([j+(t),0(t)] Jf(t)) 



o 



d ~ 
d9sm9Q{cos9)D j i [—] exp^r^ cos 6]([jf(t),0(t)] Jf(i)) 



2hc 3 J y ' Jll du. 

o 

1 

2^ 



xCM)tf#(*),o(t)] Jf(*)). (27) 



Here Tj and x(j, Z) are the spontaneous emission time and exchange integral between the single photon 
radiators respectively [13] 

the expressions in equation (|27|i and (|28|) are defined below 

A(0) = 1 + cos 2 cos 2 6>(3 cos 2 
Djl [^-] = [l + C os 2 (6) + (3cos 2 6-l)^|^] 2 ] (29) 

where cos is the scalar product between the unitary vectors along the direction of dipole momentum of 
the j (or I) atom = dj/di and the direction of the distance between the j and I atoms ny; = Tji/fji. 

2. The two-photon exchange integral between dipole forbidden transition of the radiators of D 
subsystem is described by third term in the right hand site of equation f| 13[) 

2-7T 7T OO 2lT TT OO 

t/6 y2 [ j f ao ■ a f 1.2 ji f a f ao ■ a f l.2au (neg,e Al ) 2 (n el ,eA 2 ) 2 g 2 (wi,w 2 ) 
Vj t = 7^3e y y Wisinwi y fejafci y y d9 2 sm9 2 J k 2 dk 2 2 ^ 

00 000 

t 

x / drexpH(2a;o-a;i-a;2)T + i(k 1 + k2,r j -r / )]^[A(t),0(t)]A(t-T)). (30) 







The exchange integral in Born-Marcov approximation for two-photon emission was obtained in paper [9] . 
Here we will estimate the exchange integral of expression Q30p integration firstly on the wave vectors k! 
and k 2 . Indeed considering that the amplitude q 2 (u>i, uj 2 ) is the smooth function of the variables k\ and 
k 2 in comparison with rapid oscillation functions exp[i(ki + k2, v,j — 17)] and exp[— i{2ujq — ui^ — Wfc 2 )r], 
we can approximate the amplitude q 2 (u>i,uj 2 ) according with the resonance frequencies of the atomic 
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system g 2 (wi,u;2) ~ q 2 (u>q , ujq ) . It is not difficult to observe that the maximal value of the smooth 
function under the integral is obtained for frequencies ujk 2 — ^fei ~ w o- In order to integrate this 
expression on the variables k\ and k 2 we change the variables x% = uj^ — u>i as in expression (|26|) . and 
consider that D~j~(t) = D~^{t) exp[— 2iuJot] is smooth operator. In these approximations it is obtain the 
following expression for Vj\ 

2 2 t IT TV 

Vjl = V g (2 )4°2^ 0) J dT JdOxsine! J d0 2 sme 2 D (0i)D (9 2 )exp[iLj Q r jl (cose 1 + cos9 2 )/c] 



S(t - r jX cos0i/c)<5(t - r jt cos8 2 /c) ([£>&), 0(t)]D,(t - r)) 

where expressions Dq(9\) and Dji(d/(dwo) arc defined by the expressions (|29|) . Integrating the right 
hand site of the equation pip on the solid angle and retardation, we obtain the following approximative 



expression 



where 



V£ = -^x b (j,l)([Dt(t),0(t)}Di(t)) 



1 _ 2 2 wfr&dgi / I , 1 



2rh 3 2 2 2 7rft 2 c 6 |_ W32 + wo ^31 — 

4 4dj rji dw iu)arji/c 
This exchange integral diverges, when the distance between the radiators Tji is less than the wavelength 
Ao = 2t7c/luq. In order to take in to account the value of the exchange integral for the small parameter, 
rji/Xo << 1, let us integrate the expression (|30j) taking in to account the method proposed in papers 
[9] and [12] . In this case we obtain the following expression for ReV^ 



F(j,l) = ReV^^-ld^l^-^f 



X Xj(( w fc)Xj/( w 21 - u k ) 





1 1 



W 3 1 - LO kl W 32 + U kl 



where 



Xji{u) - (1 - cos 2 £)^-f- + (1 - 3 cos 2 £) 



3. In the right hand part of this equation (|13[) the third order terms contains the resonances 
between dipole active radiators A ,B and dipole forbidden radiators D, described by the correlation 
functions (S+(t)[R+(t), 0(t)]D~(t)), ([D+(t),0(t)}Rj(t)Sf(t)) and (R+ (t)[S+ (t),0(t)}D m (t)) . Let us 
introduced the similar approximation in the chronological interaction between the atomic subsystems 
A, B, and D. The exchange integral between three atoms is represented in the similar form as in the 
expression (13T))) 



V V 2irhd s d r (2vr) 3 . /lV f 2 , /" a , ; , 
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1 

X 



/ dx\ J dx2 J dri J e?r 2 cxp[— i(u)\ — ui r )Ti — i{u>2 — u s )t2\ 



-1-10 

d d 

X Dnl ^duj'^ Dnj ^-duj'^ exp[iuj 1 r nj x 1 /c + iuj 2 r nl x 2 /c] 

([£+(<). 0{t)]Rj(t - n)Sf(t - r 2 )); (33) 

After the substitution of variables u)\ — u r — uj\ and lj-z — u s = &2, we can approximate the smooth 
amplitude uj\uj\ ^ui\Wix{w\, 0J2) with expression (oj r ) 2 ui 2 ^ui r uj s x (<*V , u> s ) . Integrals on the new variable 
u>i and u>i give the following aspect of expression (|3"3"f 

1 1 

'-2 



Vjln-as-d = 27Ti (Xj ^j^UJ 2 s (uj r ) 2 ^/uj r UJ s x{uJ s ,UJr) J dxi J dx 2 







d d 

X Dnl \.~fcJ~\ Dn i ^[^7] exp[iuj r r nj x 1 /c + iuj s r n iX2/c] 
x ([£>+(<), 0{t)]Rj(t - r njXl /c)Sf(t - r nl x 2 /c)) 

from which follows that a; l7 and X2 > . In the Born approximation the expression for V J c ln . as _ d 

Vrs-dimJJ) = -^U(j,l, m )([D+(t),0(t)}R-(t)Si(t)) 

where 

1 /2\ 2 dgdr^s^i^K^r) 3 f 1 1 



Tfcsr \>V 4:TTC 6 h 2 {(jJ 3 2+U s UJ 3 i~UJ r 

U(j,l,m) = - ( |) 2 A, [A] D nj [ Jg-] c2 [GXP ^iM ~ 1] N^/d ~ 1] , (34) 
\Z J oui s oui r uj r uj s r n jr n i 

4. Let now found the retardation in the last correlation function term of equation (| 13 p . Taking in to 
account the retardation in the rapid oscillation part of atomic operators J (t—r) = J ± (t— r) exp[±iu;(i— 
t)] where J^, u> and t are the atomic operators, transition frequencies and delay time for atomic 
subsystems A, B and D respectively, 



t t ■ (dr-,gfe 1 )(d s ,gfc 2 )(n eg ,e Al )(n ei ,e A2 )g(a; fcl ,a;A ;2 ) 

v M-d(t) = * 2^1^ 1^1^ p 

kik2 m— 1 I — 1 j=0 

t t 

x y dn exp[i(2w - uj kl - Wkjn] ^ dr 2 exp[-i(w s - )t 2 ] 


x exp[- i (k 1 ,r J - r m ) + i(k 2 ,r ; - r m )]{S+(t - r 2 )[i?+(<), 0(t)]^ m (i - n)). (35) 

Passing from the summation to integration in expression (|35l) we obtain following expression for corre- 
lation between the j, I and m atoms 
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x J dr 2 exp[-i(uj s - uj 2 )t 2 ] 
o 

d d 

xDjA- — ]Dj m [-= — 1 exp[-iu 2 r m iX2/c- iuiirj m xi/c)] 

OUJ 2 OUJ\ 

x(S+(t - T 2 )[R+(t),0(t)}D m (t - n). (36) 

Introducing the new variables u\ = u>i — w r \ u 2 = lu 2 — lu s in (|36[) . and approximating the smooth 
amplitude uj\uj 2 ^ijJiui 2 x{<^i, W2) with expression W5W^w s w r x(w s , uj r ) we get using the delta functions 

urn 

2 1 1 t t 

Vji m ;as-d(j, 1 ) = 2ii - -^Lo 2 s Lol^bj s Lo r x{u s , 0J r ) / cfei / dx 2 / (in / dr 2 



-l-io o 



x (J(t 2 - n - r J ix 2 /c)S{r 1 - r jm xi/c) 



d d 

XDjl ^dur^ Djm ^-duT^ exp[iuj s r m ix 2 /c]exp[-icj r rj m xi/c] 
x{S+(t - r 2 )[Rf(t),0(t)]D-(t - ri)), 

the value of which can be estimated observing from the arguments of ((-functions that x\ > and 
fjmXi + rjix 2 > 0. In this case, neglecting the retardation t\ and t 2 in the smooth correlation function 
we obtain 

V 3lmias . d (j,l) = -^VU,l.m)(S+(t)lR+(t),0(t)}D m (t)), 
where cooperative rate is 

1 (2\ 2 d s d r d 23 d 31 ^ujf ( 1 1 



Tsbr \3/ At7C 6 H 2 \ W32 + LO s W31 — UJ r 

the integral V(j, l.m) on the direction of the emitted photons is 

VU.Un) = D ml[ jL]D lm [JL} 

r c 2 {2exp[i(u; s r m ; - uj r r jm )/c} - 2 exp[2zo; s r m ;] - exp[~2iuj r r jm ] + 1]} ^ 

2uj s uj r ri m r m i 

, c 2 (uj s + uj r )[cxp(-iu! r rmj /c) - 1] exp(zo; s r m ;/c) - exp[-io; s r m ;/c - zuyr m j/c] 

(w s + u r )u) s w r r m jr m i 

c 2 [o; s exp[-i(cj s + w r )r mJ /c] 
+ 7 ; \ jl^^mj -Tjl). (37) 

If we will consider that X\ > 0, and x 2 > the expression for (|37p takes more simple form 

VU.Un) = (§)' 

c 2 {exp[^r mi /c] - l][exp[-^. r r jm /c] - 1] 
uj s ui r ri m r m i 

The expression for exchange integrals described in points 1. — 4. are used in the master equation (ITH1) 



References 

[1] B. Nikolaus, D. Zhang, and P. Toschek, Phys. Rev. Lett. 47, 171(1981). 



18 



[2] D.J. Gauthier, Q. Wu, S.E. Morin, and T.W. Mossberg, Phys.Rev. Lett. 68, 464 (1992); 

[3] M.Bruno, J. M. Raimond, S. Haroche, Phys. Rev. A 35, 3888 (1987); M. Bruno, J. M. Raimond, 
P. Goy, Phys. Rev. Lett. 59,1898 (1987). 

[4] G.Brcit and E.Teller, Astrophys. J . 91, 215 (1940); 

[5] M. Goeppert-Mayer, Ann Phys. (Leipzig) 9, 273 (1931). 

[6] L. Spitzer and J.L.Greenstain, Astrophys. J . 144, 215 (1940). 

[7] J. Shapiro and G. Breit, Phys Rev. 113, 179 (1959). 

[8] S.G.Alexander and P. M.Eszaros, Astrophysics J. 327, 554 (1991); 565 (1991). 

[9] N.A.Enaki, Zh. Eksp. Teor.Fiz.94, 135 (1988); N.A.Enaki, Zh. Eksp. Teor.Fiz., 98, 783 (1990); N. 
Enaki, M. Macovei, Phys. Rev. A 56, 3274 (1997). R. Marrus and R. Schmieder, Phys. Rev. A 5, 
1160 (1972). 

[10] R. Marrus and R. Schmieder, Phys. Rev. A 5, 1160 (1972). 
[11] R. H. Dicke, Phys. Rev. 93, 99 (1954). 
[12] N. Enaki, M. Macovei, Phys. Rev. A 56, 3274 (1997). 
[13] Gross M., Harochc.S. Phys.Rcpp. 93, 301-396 (1982). 

[14] Andrccv A.V., Emelyanov V.I., and Ilinskii Yu.A. Cooperative Effects in Optics (IOP Publishing, 
Bristol, 1992). 

[15] W. Heitler, The Quantum Theory of Radiation 3-ed ed., London 1954, Sec. 8] 
[16] Y. Cheng, B. Xia 2007 arXiv: 0706.0960v2. 

[17] X. Lu, J. H. Brownell, and S. R. Hartmann, Laser Physics, 5, 522 (1995). 



19 



